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Using density functional theory, we investigate the struc- 
ture of mixed 3 Hejv 3 - 4 He]v 4 droplets with an embedded im- 
purity (Xe atom or HCN molecule) which pins a quantized 
vortex line. We find that the dopant+vortex+ 4 Hejv 4 com- 
plex, which in a previous work [F. Dalfovo et al., Phys. Rev. 
Lett. 85, 1028 (2000)] was found to be energetically stable 
below a critical size N CT , is robust against the addition of 3 He. 
While 3 He atoms are distributed along the vortex line and on 
the surface of the 4 He drop, the impurity is mostly coated by 
4 He atoms. Results for N4 = 500 and a number of 3 He atoms 
ranging from to 100 are presented, and the binding energy 
of the dopant to the vortex line is determined. 
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Helium nanodroplets have recently attracted a consid- 
erable interest. A major reason is the possibility of using 
them as an inert, ultracold matrix for molecular spec- 
troscopy studies jE). They also allow one to address su- 
perfluid phenomena at a microscopic scale [|2| , and consti- 
tute an ideal testing ground for quantum many-body the- 
ories. An interesting perspective in this direction is the 
investigation of quantized vortices in finite systems. Key 
experiments have been carried out in the last two years 
on vortices in Bose-Einstein condensed gases of rubidium 
and sodium atoms confined in magnetic traps [Bj, where 
vortical states are created by acting with external per- 
turbations in different ways. These states turn out to be 
more robust than expected on the basis of qualitative ar- 
guments. In principle, analogous vortical configurations 
are also possible in superfluid 4 He droplets, where the 
external perturbation could be a moving and/or rotating 
impurity. Although the presence of vortices in superfluid 
4 He drops is not energetically favorable B, we have ar- 
gued that they can be stabilized by molecules hosted in 
the bulk of the drop || , and that their existence could be 
inferred from the changes they induce in the molecular 
spectrum ||. 

The aim of this work is to extend our previous analy- 
sis H to the case of mixed 3 He- 4 He droplets. The addi- 
tion of 3 He atoms to doped 4 He droplets has significant 
consequences in current experiments, since it lowers the 
temperature of the droplet from about 0.4 K to 0.15 K, 
and their presence can be used as another source of in- 
formation for characterizing the interaction of the dopant 



with the superfluid environment M. In this context, an 
accurate description of the first solvation layers of 4 He 
and/or 3 He around the impurity is important. The pres- 
ence of 3 He in a droplet hosting a quantized vortex would 
display several interesting features, since: i) 3 He atoms 
behave as a normal component in the superfluid, provid- 
ing a friction mechanism for the motion of vortex lines; 
ii) 3 He atoms occupy surface states, known as Andreev 
states, which are energetically favored by the lighter mass 
and the larger zero point motion of 3 He compared to 4 He; 
iii) some 3 He atoms will be attached to the vortex core, 
where they are expected to have a binding energy of the 
order of 2-3 K g] . We investigate the effects (ii) and (iii) 
by using a density functional method to calculate the 
structure and the energetics of these systems, also in- 
cluding dopant atoms and molecules. In the calculations 
we use Xe or HCN embedded in droplets with N4 = 500 
atoms of 4 He and a number of 3 He atoms N3 ranging 
from to 100. 

Our starting point is a density functional previ- 
ously developed for mixed 3 He- 4 He systems, which al- 
lows one to write the energy of the mixture as E = 
fdrH[p 3 (r),p 4 (r)}, where p 3 (r) [p 4 (r)] is the 3 He [ 4 Hc] 
particle density (see [pUlCJ] and references therein). To 
keep the already cumbersome calculations at an afford- 
able level, here we use a slightly simplified version of the 
same functional; namely, we take the core of the screened 
Lennard-Jones He-He potentials as in the original Orsay- 
Paris functional pd| , and drop the gradient-gradient term 
which appears in the Orsay-Trento functional ]15| ] for 
4 He. We have checked that these changes have negligible 
effects on the relevant results, while drastically reduce 
the numerical effort. As discussed in pJlQ|], the density 
functional contains a set of parameters which are fixed to 
reproduce static properties of pure and mixed He systems 
at zero temperature, like the equation of state, surface 
tension of the different interfaces, the osmotic pressure 
and maximum solubility of 3 He into 4 He |l2| . 

As in ||, the vortex line is included through the 
Feynman-Onsager ansatz, i.e., by adding an extra cen- 
trifugal energy associated with the velocity field of 4 He, 
which is singular on the vortex axis, thus forcing its den- 
sity to vanish. For doped droplets, one has to include 
the helium-impurity interaction, which acts as an exter- 
nal potential in which the helium density adjusts to min- 



imize the energy. The potential for Xe has been taken 
from Q, and that of HCN from |0|. Combining all 
terms, the total energy can be written in the form 



E = Jdr jw[p 3 (r),p4(r)] + 
+ y / (r)[p 3 (r)+p 4 (r)]} , 
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where Vi is the helium-impurity potential and r± is the 
distance from the vortex axis. The energy minimization 
is performed in axial symmetry by mapping the densities 
on a spatial mesh, putting the vortex line along the z-axis 
and the dopant in the center, at r = 0. 

First we consider two rather simple configurations: i) 
a mixed droplet with a dopant, but no vortex; ii) a mixed 
droplet with a vortex, but no dopant. In the former case 
we just re-obtain the results of d|H]], namely that the 
amount of 3 He atoms in the bulk of the drop is negligible, 
and the dopant is coated by 4 He. The structure of the 
drop is 'onion like', with 3 He distributed in the outer shell 
and He inside, surrounding the embedded impurity. In 
contrast, in the latter case 3 He atoms can funnel through 
the vortex dimple created at the 3 He- 4 He interface and 
eventually the vortex core is filled with 3 He. 

An example of a mixed droplet with a vortex is shown 
in Fig. [j]. The 4 He and 3 He density profiles in the radial 
direction, at z = 0, arc shown in Fig. Bfor different values 
of A3. It can be seen that 3 He is filling the vortex core 
even for A3 = 20, and that the dependence of the central 
density on A3 is weak, as expected for a close-packed lin- 
ear chain of atoms. The remaining 3 He atoms occupy the 
available surface states, whereas 4 He stays in the bulk. A 
comparison with the case of pure 4 He droplets (solid line 
in Fig. pf) shows that 3 He atoms push the superfluid 4 He 
component away from the vortex axis, thus lowering the 
kinetic energy associated with the vortex flow |16],[l7| ■ 

Placing a dopant in the center of the droplet signif- 
icantly distorts the 3 He and 4 He densities. In Fig. || 
we show the same densities as in Fig. ll but with an 
embedded HCN molecule. The strong helium-impurity 
attraction, which is the same for 3 He and 4 He, favors 
the formation of a layered structure of 4 He atoms near 
the molecule, since they have a smaller zero point mo- 
tion than 3 He atoms and can be localized more easily 
in the local minima of the potential. In the vicinity of 
the dopant, 3 He atoms only remain at the pinning points 
on the vortex axis, where 4 He is excluded from by the 
high kinetic energy of the vortex flow. The net effect is 
that the structure of the complex dopant+vortex near 
the dopant is very similar to that of pure He droplets. 
It is worth to see that the dopant also produces a mod- 
ulation of the 3 He density along the vortex line. 

Let us denote with subscripts X and V the energies, 
E, of mixed droplets doped with an impurity X and/or 
containing a vortex line. The energetics of these droplets 
can be studied by introducing, for a fixed value of A4, 



the following functions of A3 [pj : 

AE V (N 3 ) = E V {N 3 ) - E{N 3 ) (2) 

AE*(N 3 ) = E X+V {N 3 ) - E X {N 3 ) (3) 

S X (N 3 ) = E X {N 3 ) - E{N 3 ) (4) 

S X +v{N 3 ) = E X+V (N 3 ) - E(N 3 ) (5) 

S X (N 3 ) = AE$(N 3 ) - AE Y (N 3 ) . (6) 

The quantities AEy and A£y correspond to the en- 
ergy associated with the vortex flow in a droplet with- 
out dopant and with dopant X , respectively. The quan- 
tity Sx is the solvation energy of impurity X in the 
mixed cluster, while Sx+v i s the solvation energy of the 
dopant+vortex complex. When the quantity S x is nega- 
tive, its absolute value represents the binding energy of 
the dopant to the vortex in the mixed cluster. 

In Ref. M we have studied the above energies in the 
A 3 = case, finding that S x+ y(N 3 = 0) is negative and 
hence the 4 Hc+X+vortex complex is stable for values of 
A4 smaller than a critical number, N cr , of the order of 
8000 for both Xe and HCN. The effect of a non-zero A 3 
value comes from a delicate interplay between different 
energy contributions, which are sensitive to the distribu- 
tion of 3 He atoms in the vortex and near the dopant. The 
crucial question is whether 3 He may depin the impurity 
from the vortex core, i.e., S X (N 3 ) becomes positive for a 
certain N 3 . 

The relevant energies are plotted in Fig. H. The top 
panel shows how the vortex energy decreases with A3 in 
droplets without dopant, with a Xe atom, and with a 
HCN molecule. This behavior is consistent with a reduc- 
tion of the kinetic energy of the vortex flow when normal 
3 He atoms displace superfluid 4 He atoms away from the 
vortex core. The two middle panels show the solvation 
energy of the impurity and of the dopant+vortex com- 
plex; both are negative for these droplets. The main re- 
sult of this analysis is shown in the bottom panel, where 
one may see the effect of 3 He on the binding energy of 
the dopant to the vortex. The binding energy \Sx\ of the 
dopant decreases when N 3 increases, but the dopant is 
still pinned to the vortex. The initial slope of the curve is 
steeper than for large A3 . This is consistent with the first 
3 He atoms occupying states along the vortex line close to 
the dopant, thus affecting 5x in a more significant way. 
In the intermediate region for A3 = 20-40, the binding 
energy exhibits a plateau. In this range, the presence of 
the dopant is expected to affect the distribution of 3 He 
atoms both in the vortex line and at the surface. One 
can see in Fig. that in the same range of A3 values 
the mixed droplet starts to develop a 3 He 'skin', i.e., an 
outer shell where the 3 He density is larger than the 4 He 
density. Whereas its effect is imperceptible at the scale 
of the energies defined in Eqs. g-|5|, it shows up in Sx, 
which is around two orders of magnitude smaller Jig ] . 

We conclude that the dopant+vortex+ 4 Hejv 4 complex 
is robust against the addition of moderate amounts of 



3 He atoms. This may offer some experimental advan- 
tages. On the one hand, mixed droplets reach lower tem- 
peratures than pure 4 He droplets; on the other hand, 
adding a variable amount of the normal component will 
result in a variable damping for the vortex motion, with- 
out loosing the characteristics that make 4 He drops ap- 
pealing for molecular spectroscopy, since the dopant en- 
vironment essentially consists of 4 He atoms as in pure 
drops. 
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means that it is low enough to render thermal excitation 
effects negligible in both 3 He and 4 He, but it is above the 
critical temperature for 3 He superfluidity. 
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FIG. 1. Density distributions of 3 He (top) and 4 He (bot- 
tom) in the xz plane for the 4 He5oo+ 3 Heioo droplet hosting a 
vortex line along the z axis. Lengths are in units of A. Darker 
regions are high density regions. 



FIG. 2. Density profiles of 4 He and 3 He in the radial 
direction, at Z — 0, for droplets with a vortex line along the z 
axis and with N A = 500 and N 3 = 0, 20, 50 and 100. The 3 He 
density profiles appear in two disconnected parts separated 
by the corresponding 4 He density profile. For A3 = 100, the 
dashed lines correspond to a cut at z = of the densities in 
Fig. hi. Lengths are in units of A. 

FIG. 3. Same as in Fig. hj but with a dopant HCN molecule 
in the center of the droplet. 



FIG. 4. From top to bottom panel: Vortex energy of the 
4 He5oo+ 3 Hejv 3 ; solvation energy of Xe and HCN dopants; sol- 
vation energy of the dopant+vortex complex; binding energy 
<5x|. The triangles represent results for Xe, the squares for 
HCN, and the circles in the top panel are the results for un- 
doped droplets. The energies are in units of K, and the lines 
have been drawn to guide the eye. 
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